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The diffraction effects on gravitational waves propagating through a stellar cluster are 
analyzed in the relevant approximation of Fresnel diffraction limit. We find that a grav- 
itational wave scintillation effect - similar to the radio source scintillation effect - comes 
out naturally, implying that the gravitational wave intensity changes in a characteristic 
way as the observer moves. 

Keywords: gravitational waves: general 

1. Introduction 

General Relativity predicts the existence of gravitational waves that propagate in 
the vacuum with the speed of light. Analogously to electromagnetism, it i s expe cted 
that gravitational waves might also be scattered by inter vening particles I ' " I and 
micro-lensed by compact objects along the line of sight l 4 l v l v l Moreover, gravita- 
tional waves will be diffracted as they pass through a distribution of objects, each of 
which acts as an obstacle (or slit). Of course, these effects are expected to become 
apparent only if the linear sizes of the intervening particles are comparable to the 
gravitational radiation wavelength in the range 10 7 — 10 15 cm. We consider a stel- 
lar cluster as the diffracting object and concentrate in particular on that possibly 
hosted in the center of the Milky Way and surrounding the central black hole Sgr 
A*. 

There are two types of diffraction phenomena discussed: the Fraunhoffer t3 and 
Fresnel approximation. In our case only the second is relevant since the first approxi- 
mation involves infinite source-cluster and cluster-observer distances. We determine 
the conditions under which the diffraction effects may be measurable. 

It turns out that a gravitational wave scintillation effect, analogous to the well 
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known scintillation of radio sources, comes out naturally implying that the grav- 
itational wave intensity changes in a characteristic way with the motion of the 
observer. 

The paper is structured as follows: in Section 2 we describe the parameters of 
the cluster mass density profile used to study the Fresnel diffraction. In Section 3 
is given the mathematical tools of the Fresnel theory, in the particular geometry 
described there. In Section 4 we present the main result of this paper: namely, 
the scintillation pattern of the gravitational waves as they pass through the stellar 
cluster at the galactic center. Finally, in Section 5 we present some conclusions. 



2. The Mass Distribution Model 

The ESO and Keck teamsEHhave continuously monitored the region of the galactic 
center for several years so that a cluster of stars (at a distance < 1") has been 
identified. In particular, Ghez et al. have reported on the observations of the 
main sequence S2 star (with mass M$2 — 15 Mg) orbiting the central black hole 
with a Keplerian period of ~ 15 yrs. This has allowed 121 the mass contained within 
R ~ 4.87 x 10~ 3 pc to be constrained to M(R) ~ 3.67 x 10 6 M Q . The most plausible 
model is that this mass is in the form of a super-massive black hole. 

However, there is a possibility that a small fraction of the inferred central mass 
is in the form of a diffuse stellar cluster which surrounds the central black hole. In 
this case, according to the dynamical observations towards the galactic center, we 
require that the total mass satisfies the condition M(R) — Mbh + Mcl(R), Mbh 
and Mcl(R) being the black hole and the cluster mass within R, respectively. 

Of course, the cluster mass and density distribution are still unknown so that, as 
a toy model, we can assume that the stellar component follows a Plummer density 
profile given by 



Pcl{t) = pof(r) , with f(r) 



1 



-5/2 



(1) 



po = rRcL . „.„ : » ( 2 ) 



where the cluster central density po is given by 

M CL 
J" CL 4nr 2 f{r) dr 

Rcl being the cluster radius. Useful information is provided by the black hole mass 
fraction, X BH = M BH /M and its complement, X C l = 1 — ^bh- As one can see, 
the requirement given in eq. J5J implies that M(r) — > Mbh for r — > 0. Hence, the 
total mass density profile p(r) is given by 

p(r) =X BH MS^(T')+pof(r) (3) 

and the mass contained within r is 

M(r) = X BH M + I 4irr' 2 p f(r') dr' . (4) 
Jo 
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Considerations on both the dynamics and evolution of the stellar cluster surround- 
ing the central black hole allow us to constrain the possible range of values El f or 
^bh to be above about 0.9. In the following for definiteness we will assume that 
\ BH ~ 0.99. 

Assuming that all the cluster components have the same mass (~ 1 M©), 
we can estimate the average distance between two neighboring stars as (d(r)) ~ 
[M Mr)]V3. 



3. The Fresnel Diffraction Approximation 

Consider a distant source emitting gravitational radiation. Assume that there is a 
distribution of objects between the source and the observer, each of which acts as 
an obstacle (or slit). Let (Xo,Yo), (X±,Yi) and (J^,^) be three reference frames 
(with mutually parallel axes) in the planes of the observer, diffraction and source, 
respectively chosen such that the origins are collinear. 

Using the Fresnel approximation for interaction of ele ctromagnetic waves with 
dust particles distributed along the line of sight givesLJ-U^the diffraction amplitude 
Aq on the plane of the observer (at distance zq from the diffraction plane) is 

p iKZ r poo 

A (x ,yo) = ^J J A , 1 (x 1 ,y 1 )e^ [{x °- Xl) +{v °-^ ] ] d Xl dy u (5) 

where 

A p iwi2 

A[(x uyi ) = ^ e **(«i.vO (6) 

is the amplitude just past the diffraction plane. In eq. 10 A2 is the source amplitude, 
r\i ~ z\ is the distance of the source from a generic point on the diffraction plane, 
and z\ is the distance of the source from the center of the diffraction plane. It is 
worth noting that when the wave moves through the diffraction plane, a change 
of phase does occur due to the optical path difference 5(xi,y±) which, in turn, 
depends on the coordinates of each obstacle. For incoming plane waves generated 
by a point-like source, the factor A — A 2 e 1Kri2 / ' Z\ is a real constant so that the 
previous integral becomes 



Ae l 



■ jxQ-xi) +(yp-yi)- J 



Aoixcyv)^— W J J e iK5 ^^e ^ dx x dy u (7) 

where the Fresnel radius is defined as 

It is assumed l ^ l ^ l that the screen is a step of optical path 6 parallel to the 
Y\ axis, described by a Heaviside distribution in the {Xi,Y\) plane 8{x\,yi) — 
5 x H(x\). This case is realistic since, at the Fresnel scale, the edge of a particle 
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distribution can be considered as a straight line that divides the plane into two 
regions. Under these assumption, the integral in eq. Q turns out to be 



A (x ,y ) = Ae 1 ' 



1 + + S(X ) + C(X ) - i(S(X ) - C(X ))} 



(9) 



where S(X) and C(X) are Fresnel integral and X$ = xq/(^JitRf) is the reduced 
variable corresponding to xq. Thus, the intensity on the plane of the observer is 

Io(xo,Vo) = A Q (xo,yo)Ao(x ,yo) = A 2 i (X ) , (10) 



io(Xo) being the relative intensity given by 



i (X ) = 1 - (S(X ) - C(X ))sin( K 5) + 



S(X ) 2 +C(X ) 2 -± 



[I-cos(kS)] .(11) 



In the next Section, we apply the previous formalism to the case of monochro- 
matic gravitational radiation (with wavelength Xgw) passing through the stellar- 
cluster possibly surrounding the galactic center black hole. Taking the optical path 
6 to be of the same order of magnitude as the average distance (d(r)) between the 
star in the cluster, we give the condition under which Fresnel diffraction occurs. 
Moreover, we shall find that a characteristic scintillation pattern appears in the 
plane of the observer. 

4. Gravitational Wave Scintillation 

The Fresnel diffraction limit entails that D>Rp, where D is the circular hole 
diameter and Rf is the Fresnel radius. In order for there to be diffraction by the 
stellar cluster, the optical path 6 should be of the same order of D that is ~ (d), 
i.e. the average distance between neighboring stars. Therefore, we get 

V Zn 

implying 

X GW < 2 -^. (13) 

zo 

Thus, in order to have diffractive effects by the star cluster in the Fresnel limit, we 
need gravitational radiation with wavelength less than the mean spacing squared 
divided by the distance of the observer from the plane of diffraction that we take 
to be zq ~ 8 kpc. 

Next we consider a gravitational wave source behind the stellar cluster at the 
galactic center. The source might be, for example, an extragalactic system of in- 
spiral binary black holes or a single spinning pulsar or a binary system of neutron 
stars at the outer edge of the Galaxy, which are assumed to be detectable by a 
gravitational wave telescope (like VIRGO, LIGO, LISA or ASTROD) with an in- 
tegration time T (for details see e.g. Thorne^j). In the following, we shall mainly 
concentrate on diffraction by the star cluster at the galactic center and, due to 
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the expected average distance between neighboring stars, the most relevant gravi- 
tational radiation sources are pulsars lying behind that cluster. Therefore, we will 
consider in the following analysis the VIRGO interferometer. 

In Fig. ^ we present the scintillation pattern obtained from eq. Ill If in the 
plane of the observer for an incoming monochromatic gravitational wave whose 
wavelength is Xqw — 1-5 X 10 7 cm. The diffracting stellar cluster has parameters 
^bh = 0.99, po = 1.15 x 10~ n g cm~ 3 and r c — 5.8 mpc, corresponding to an 
average stellar distance of (d) ~ 10 15 cm evaluated at a distance of r ~ r c from the 
cluster center. The position Xq = corresponds to the case of perfect alignment 
among the source, the cluster center and the observer. Of course, as one can see 
from Fig. El as the observer moves along the direction Xq, it passes through a series 
of maxima and minima of the gravitational wave intensity. 



j , 

-4-0 4 
Xo 

Fig. 1. The gravitational wave pattern in the plane of the observed is shown. Here, the diffracting 
stellar cluster has parameters Xbh = 0.99, po = 1.15 X 10 -11 g cm -3 and r c = 5.8 mpc, 
corresponding to an average stellar distance of (d) ~ 10 15 cm evaluated at a distance of r ~ r c 
from the cluster center. 

Analogously to Fig. ^ the gravitational wave pattern obtained with (d) ~ 10 14 
cm, evaluated at the distance r ~ 0.1r c , is reported in Fig. [3 

As one can see by comparing Fig. [21 with Fig. ^ even though the frequency of 
the curves remains the same, the height of the curves decreases and the principal 
maximum shifts from «o — 1-8 to iq ~ 1.4. If one considers a different gravitational 
radiation wavelength, which is to say a different wave source, we obtain analogous 
scintillation patterns to those presented in Fig. ^ (for (d) ~ 10 14 cm) or Fig. [21 (for 
(d) ~ 10 15 cm). 

Obviously, the positions of the maxima do not change in the dimensionless vari- 
able Xq. However, different gravitational radiation diffraction patterns (for different 
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Fig. 2. The same as in Fig. 1 but considering, as diffracting objects, the stars with an average 
distance (d) ~ 10 15 cm evaluated at r ~ 0.1r c from the cluster center. 



incoming wavelengths) are obtained in the variable xq. As is clear, each figure has 
been normalized to the incoming gravitational radiation intensity (io = 1)- As such, 
as the observer moves, the scintillation effect implies a change in the wave intensity 
that, in principle, might be used to infer information about the source and the 
cluster parameters (in particular on (d)). 

As an example, consider a gravitational wave (with wavelength \qw — 1.5 X 10 7 
cm) interacting with the stellar cluster at the center of the Galaxy. By inspecting 
Fig. 1, it is clear that the distance between the principal maximum and the next 
minimum is AATo ~ 1.5. Hence, if the observer moves around the center of the 
Milky Way with the typical orbital speed Vq = 220 km s^ 1 , the time-scale AT ~ 1 
year is obtained. For such a cluster we do not consider binary sources because of 
the Fresnel diffraction limits |JT3}. Obviously, if we take into account other types 
of clusters (with different physical parameters) it is also possible to get diffractive 
effects with different gravitational wave sources, like black hole binaries, etc. 

In order for the diffraction effect to be measurable, it is clear that the integration 
time T, required by a gravitational wave interferometer to detect a signal over the 
typical instrument noise, has to be less than the traveling time, Tq, for Earth to 
move from a maximum to the next minimum of the scintillation pattern. Here, an 
optimal detection sensitivity or, roughly speaking a sig nal-to-noise ratio ~ 1, is 
assumed. Consequently, the integration time becomes 



T 



~Hf)~ 



(14) 
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where h(f) is the experimental sensitivity curve (in unit of y/ Hz) and ho is the ex- 
pected (dimensionless) gravitational wave amplitude at the characteristic frequency 
/■ 

In Tab. Q we give the relevant parameters for three pulsars and the resulting 
values for T© and T to be evaluated by taking into account the VIRGO sensitivity 
threshold a . As one can see, for the Crab pulsar the scintillation effect is clearly 
detectable, for PSR 0021-72C the effect is marginally detectable while for PSR 
1604-00 it is far from detectability, at least for the VIRGO interferometer. 

Table 1. Here the physical parameters for some known pulsars together with the characteristic 
frequency / at which the objects mostly emit are reported. The expected gravitational wave am- 
plitude ho and the experimental sensitivity curve h(f) (taken from VIRGO sensitivity threshold) 
at that frequency are also reported. The integration time T and the traveling time Tq are given 
with the Jninimum average distance (d) between two neighboring stars (see eg. 1131 . 



Pulsar 


Crab 


PSR0021-72C (47TUC) 


PSR1604-00 


P (ms) 


33.2 


5.757 


421.8 


P xlCT 15 (s/s) 


422.6 


0.04 


0.30 


/ (Hz) 


60.423 


347.403 


4.742 


ho 


3.59 x 10" 2b 


8.36 x 10~ 27 


2.68 x 10" 27 


Hf) 


7 x 10~ 23 


5 x IO" 23 


2 x 10~ 20 


Tq 


~ 3.75 y 


~ 1.563 y 


~ 13.38 y 


T 


~ 10.6 h 


~ 1.1 y 


~ 2 X 10 6 y 


(d) (cm) 


1.4 x 10 ib 


5.68 x 10 ib 


5.98 x 10 14 



In Tab. we give the same quantities as in Tab. ^ but for three different kinds 
of pulsar parameters: a typical radio pulsar, a typical new-born pulsar and a mil- 
lisecond pulsar. As one can see, only for new-born pulsars, even with rather small 
P values, the scintillation effect would be detectable by the VIRGO detector. 

Table 2. This tabli is analogous to Tab.^but for three different classes of pulsar parameters. 



P(s) 


0.1 


0.01 


0.001 


P (s/s) 


io- 15 


io- 15 


10 -w 


/ (Hz) 


20 


200 


~ 2000 


ho 


~ 10" 2b 


~ 3 x IO" 26 


~ IO" 27 


Hf) 


5 x 10" 22 


4 x 10 -23 




Tq 


~6.5y 


~2y 


~ 0.92 y 


T 


~79y 


~ 18.5 d 


~ 300 y 


(d) 


2.42 X 10 15 


7.66 x 10 14 


2.42 X 10 14 



a See VIR-NOT-PER-1390-51 available online at |http: //www.virgo.infn.it/senscurve| 
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5. Concluding Remarks 

The General Theory of Relativity predicts the existence of gravitational waves. In 
accordance with what is currently observed for electromagnetic waves interacting 
with obstacles (or slits), we expect that also gravitational radiation may undergo 
diffraction. As usual, it is possible to have diffraction effects only if the incident 
radiation wavelength and the intervening obstacle size are of the same order of 
magnitude. Hence, due to the typical large wavelength of the considered radiation 
(i.e. 10 7 — 10 15 cm for spinning neutron stars), only systems of objects on galactic 
scale may act as diffractive objects. We have considered the diffraction of grav- 
itational waves by a stellar cluster possibly hosted at the center of our Galaxy 
and found, under the conditions for which the Fresnel diffraction occurs, that a 
scintillation of the gravitational wave signal comes out naturally. 

As we have shown, this effect, analogous to the well known scintillation of radio 
waves due to a medium along the line of sight of a distant source, produces a 
characteristic pattern on the plane of the observer. We have also shown that the 
typical crossing time Tq (i.e. the time required by the observer to move across a 
minimum to the next maximum of the gravitational wave scintillation pattern) is of 
the order of a few years. For a particularly strong gravitational wave signal and in 
particular regions of the frequency band, the typical integration time of instruments 
like VIRGO and/or LIGO (which have a rather similar sensitivity threshold curve) 
may be as short as few hours thus allowing a good sampling of the scintillation 
pattern. 

We emphasize that, as observed in Section 4, new-born pulsars, i.e. pulsars 
that were born within the last ~ 10 4 years are the best candidates to look for the 
gravitational wave scintillation effect by the star cluster lying at the galactic center 
with the VIRGO instrument. 

It is worth noting that the next generation of gravitational wave interferometers 
(such as Advanced VIRGO and Advanced LIGO) will have a sensitivity greater 
than that allowed by the present instrument configurations. For example, Advanced 
LIGO b will have more than an order of magnitude greater sensitivity than initial 
LIGO implying a reduction by at least two orders of magnitude of the integration 
time. If this is the case, there will be a real possibility of using the gravitational 
radiation signals to probe deep into our Universe and search for effects like the 
scintillation of gravitational waves. Detecting this effect will give an independent 
way of estimating the main parameters of the intervening stellar cluster like the 
average stellar distance (d) and the cluster central mass density. The microlensing 
effect on the gravitational waves due to the black hole at the galactic center has 
been investigated elsewhere^. 



b See the web page at http://www.ligo.caltech.edu/advLIGO/ 
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